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Summary of Technical Progress 


Objectives ■ The overall objective of this research is to develop models to predict drop-size- 
distribution evolutions due to droplet collisions and coalescence during ^processing within 
the miscibility gap of bimetallic liquid-phase-miscibility-gap materials. The individual and 
collective action of gravitational and nongravitational mechanisms on the relative motion 
and coalescence of drops are considered. 


Research Task Description, : When bimetallic liquid-phase-miscibility-gap materials, which 
are thought to have a variety of desirable properties, are cooled through the miscibility gap, 
droplets rich in one of the metals form in the liquid matrix rich in the other metal. Droplet 
coalescence and phase segregation then occur due to buoyancy and to thermocapillary and 
other nongravitational mechanisms. In order to gain a predictive understanding of these 
phenomena, population dynamics models are used to follow drop-size distribution evolutions 
in time as the droplets grow due to collisions and coalescence. Continuous drop size dis- 
tributions are discretized into a large number of categories. Drops of a given mass are 
destroyed by coalescing with other drops and are formed by the coalescence of smaller drops . 
The population dynamics model tracks the formation and destruction of drops in each size 
category. 


The relative motion of drops which gives rise to their collision and coalescence is considered 
to occur by gravity sedimentation, Marangoni migration, Brownian motion, and bulk. flow. 
These collision mechanisms are considered either individually or collectively in the pop- 
ulation dynamics models. In general, different mechanisms dominate for different pro- 
cessing conditions, materials properties, and drop-size ranges. 


The collision kernals appearing in the population dynamics models require expressions for 
the collision rate between drops of two different sizes. Classical expressions attributed to 
Smoluchowski are improved to include attractive, repulsive, and hydrodynamic interactions 
between drops. In particular, trajectory calculations are used to predict collision efficiencies, 
which represent the ratio of the collision rate with these interactions to the Smoluchowski 
collision rate without these interactions, as functions of the size ratio, viscosity ratio, and 
other relevant dimensionless parameters. 


Accomplishments: During the three years of NASA support, the following progress was 
made: 

1 A computer program has been completed for solving the population dynamics model 
to follow droplet size evolutions with time in homogeneous dispersions due to col- 
lisions arising from gravity sedimentation, Marangoni migration, and/or Brownian 
motion. Some of the key results are that a bimodal initial distribution will exhibit 
much more rapid coalescence due to gravity sedimentation or Marangoni migration 
than will a unimodal initial distribution (Figure 1), a unimodal initial distribution 
will evolve into a bimodal distribution and then into a shifted and broadened uni- 
modal distribution (Figure 2), and that coalescence may be greatly reduced by 
antiparallel alignment of the gravity vector and the temperature gradient (Figure 
3). 

2. Collision efficiencies for Brownian motion (Figure 4) and gravity sedimentation 
(Figure 5) for drops having a range of viscosity and radius ratios have been computed 
both in the presence and absence of attractive forces. A key result is that, in contrast 
to rigid particles, spherical liquid drops have nonzero collision efficiencies in the 
absence of attractive forces. 

3. Theoretical work has been completed to predict collision efficiencies for Marangoni 
migration of two drops of different sizes in a temperature gradient Example results 
are shown in Figure 6. Note that the collision efficiencies for thermocapillary motion 
are larger than those for gravity motion, because of the more rapid decay of the 
velocity fields created by drops undergoing thermocapillary motion, provided that 


2 



the thermal conductivity of the drops is not too large. Surprisingly, however, the 
relative mobility function along the line-of-centers can become negative for large 
thermal conductivities (Figure 7), indicating that the smaller drop moves away from 
the larger drop and that there is a region of closed trajectories which prevents 
coalescence in the absence of attractive forces (Figure 8). 


4. In related experimental work, drop size distributions of butyl benzoate in water are 
followed with time by holography as coalescence occurs due to gravity motion. We 
are also examining the interaction of two drops with video microscopy. An example 
of coalescence is shown in Figure 9, and measured relative trajectories of two drops 
are shown in Figure 10 to be in good agreement with theory. 

5 Theoretical work on combined mechanisms for drop coalescence has been focused 
on buoyancy and Marangoni motion with the gravity vector aligned parallel or 
antiparallel to the temperature gradient. Our results show that there is a finite 
region in parameter space for which no collisions occur (Figure 11). This is because 
the interaction of two drops due to thermocapillary migration decays more rapidly 
(as 1/r 3 ) with the separation distance (r) than does the interaction of two drops in 
gravity motion (which decays as 1/r). As a result, the large drop which moves toward 
I small drop below it due to gravity may reach a separation distance where the 
initially weaker thermocapillary effects just balance the buoyancy-dnyen relative 
motion, and so the separation distance would then no longer decrease (Figure 1Z). 
With appropriate antiparallel alignment of buoyancy and Marangoni velocities, a 
polydisperse suspension of small drops is predicted to evolve into a monodisperse 
suspension of larger drops (Figure 13). 

6. Time scales have been determined for phase segregation and drop coalescence, 
together with criteria for predicting whether or not significant coalescence wi 11 occur 
prior to phase segregation. A computer model has been developed to solve population 
balances for drop size distributions which vary in space due to buoyancy or ther- 
mocapillary motion. Results show that the phase segregation rate initially increases 
due to coalescence and subsequently decreases as the larger drops migrate out ot 
suspension (Figure 14). 

7. Physical data on a variety of bimetallic and transparent immiscibles have been 
collected, and dimensionless parameters representing various effects have been 
tabulated as functions of drop size. Data on composite Hamaker constants (for van 
der Waals attractions) and temperature coefficients of interfacial tension (for 
thermocapillary migration) are sparse, and order-of-magnitude estimates are typ- 
ically used. Calculations of collision efficiencies have been made for bismuth drops 
in a zinc melt, for ethyl salicytate drops in diethylene glycol, and for lead drops in 
an aluminum melt. 
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Figure la. Time evolution ofthedrop size 
distribution for a bimodal initial distri- 
bution of bismuth drops undergoing 
Marangoni migration in molten zinc. 



Figure 2. Time evolution of the drop size 
distribution for a unimodal initial distri- 
bution of water drops subject to gravity 
sedimentation in air. 


Figure lb. Time evolution of the drop size 
distribution for a unimodal initial distri- 
bution of bismuth drops undergoing 
Marangoni migration in molten zinc. 
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Figure 3. Mass average drop radius at t»46 
seconds due to combined gravity and 
Marangoni motion of bismuth drops in 
zinc as a function of gravity strength 
relative to normal for parallel (£>0) and 
antiparallel (£ <0) orientation of gravity 
vector and temperature gradient. The 
initial mass average radius is r g *2 pm. 
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Figure 4. Collision efficiencies for viscous drops of various viscosity ratios ((I) and 
size ratios (X) undergoing Brownian motion in the absence of attractive 
forces. 



Figure 5. Collision efficiencies for viscous drops of various size ratios (X) and viscosity 
ratios (jl) undergoing gravity sedimentation in the absence of attractive forces. 
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Figure 


a Collision efficiencies as a function of the drop size ratio (X) for drops of various 

6 ' SSff«SS5S&“ »o. <i) and vi.co.ity ratio, ft.) underaouw thermoc 
p illar y migration in the absence of attractive forces. 
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Relative mobility function versus the dimensionless separation distance for 
two unequal drops (X « 0JS) in relative motion along their Une of centers due 
to thermocapillary migration, for various thermal conductivity ratios (k) and 

viscosity ratios (it). 
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Figure 9. Time sequence of approach and coalescence due to gravity of two castor oil 
drops in a silicon fluid. 



Fiaure 10 Dimensionless relative trajectories of two different sized castor oil drops 
figure iu. due to grotty m a silicon fluid. The symbols are experimental data 

and the dashed lines are theory. 







The collision efficiency versus the ratio of buoyancy and th nTfnr vaHoS 
relative velocities for two widely spaced drops of size ratio X - 0.5 for various 
viscosity (ji) and thermal conductivity ratios (*). 



dimensionless separation distance 

Figure 12. The dimensionless gravitational (L), thermocafil^ ^N^nd combined 
relative velocities versus separation distance for two drops with X - 0.5, n - 1. 

k. = 1, and N c = -2.5. 
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Fimire 13. Evolution of the size distribution of lead drops in an aluminum for 

lgU an tiparallel alignment of the gravity vector (9.8 m/s*) and temperature gra 

dient (ltf K/m). 



Fimire 14. The rate of phase separation versus time for a dispersion of 

^ undergoing simultaneous buoyancy motion and coalescence in a zmc melt m 

containers of increasing height (left to right), the dashed curves are the cor- 
responding results in the absence of coalescence. 
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